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ABSTRACT 



Extensions of Einstein gravity which allow the gravitational constant G to change 
with time as the universe evolves may provide a resolution to the horizon problem without 
invoking a period of vacuum domination and without the subsequent entropy violation. 
In a cosmology for which the gravitational constant is not in fact constant, the universe 
may be older at a given temperature than in a standard Hot Big Bang universe; thus, 
larger regions of space could have come into causal contact at that temperature. This 
opens the possibility that large regions became smooth at some high temperature without 
violating causality. The extra aging of the universe can be accomplished by an early 
period with a large Planck mass, a period we call the MAD era (Modified Aging era or 



the Massively Aged and Detained era). We discuss in this paper theories of gravity in 
which the gravitational constant is replaced with a function of a scalar field. However, 
this resolution to the smoothness problem can more generally be a feature of any physics 
which allows the Planck mass to vary with time. In this paper, we examine scalar theories 
of gravity without a potential for the scalar field. We first consider the original Brans- 
Dicke proposal and then address more general scalar theories. 

Solutions to the equations of motion for Brans-Dicke gravity during the radiation 
dominated era are presented. In particular, we study the evolution of the Brans-Dicke 
field $ which determines the Planck mass at any given time, = mpi{t)'^. We find that, 
regardless of initial conditions for the Planck mass, it evolves towards an asymptotic value 
rhpi = The same asymptotic behavior is found in more general scalar theories. For 

both a Brans-Dicke cosmology and a more general scalar theory, the smoothness of the 
universe can be explained if the Planck mass is large at some high temperature Tc prior 
to matter-radiation equality: specifically, if rhpi/Mo > Tc/Tg, where Mg = 10^^ GeV is 
the Planck mass today and Tg is the temperature of the cosmic background radiation 
today. In a pure Brans-Dicke cosmology, an additional mechanism (e.g. a potential for 
$) is required to drive the Planck mass to the value Mg by today. In a more general 
scalar theory of gravity with variable Brans-Dicke parameter, the suggestion is made that 
the Planck mass may approach the value Mg more rapidly during the matter dominated 
era than in a Brans-Dicke cosmology. 
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I) Introduction 

The standard Hot Big Bang model of the early universe is unable to explain the 
smoothness or flatness of the observed universe. In the standard cosmology, our present 
horizon volume would envelop many regions which were causally disconnected at earlier 
times. Consequently, the homogeneity and isotropy of the observed universe is a mystery. 
Regions which could not have been in causal contact at earlier times seem nonetheless to 
be identical in temperature and other properties, as the isotropy of the cosmic background 
radiation attests. 

The inflationary model proposed by Guth^ addresses the horizon and flatness prob- 
lems, as well as the monopole problem (if it exists). As a general class of early universe 
models, inflation suggests that our universe passes through an era of false vacuum domi- 
nation during which the scale factor grows exponentially (or at least super luminally) . The 
superluminal growth of the scale factor inflates a region which was initially subhorizon- 
sized and therefore in causal contact. If the scale factor grows sufficiently, our observable 
universe fits inside one of these blown up causally connected volumes. During infiation, 
the temperature of the universe plunges, T oc R~^, where R is the scale factor. There- 
fore, the next crucial ingredient for a successful inflationary model is a period of entropy 
violation which reheats the universe to some high temperature. 

In this paper, we propose that a cosmology with a variable Planck mass can resolve 
the horizon problem without a period of vacuum domination. Further, entropy is always 
conserved. [In another paper, we illustrate how our model can resolve the monopole prob- 
lem as well; the flatness problem may be alleviated by our model but is quite complicated 
and will be studied elsewhere.] We have considered (a) (in Section III) the Brans-Dicke 
proposal to replace the constant Planck mass with a scalar fleld, rupi oc ^, and (b) (in 
Section IV) more general scalar theories where the Planck mass could be an arbitrary 
function of a scalar field '0. In both cases, the energy density of the universe begins radi- 



ation dominated and then goes over to a period of matter domination as in the standard 
cosmology. 

We derive below the analytic solutions to the cosmological equations of motion for 
these alternate theories of gravity when the energy density in ordinary matter is radiation 
dominated. We find the scale factor, the temperature, the Hubble constant, and the 
horizon radius in terms of the variable Planck mass. We also find the time evolution of 
mpiit) for early and late times during the radiation dominated era. 

The initial time derivative of the Planck mass, drupi/dt, may be positive, negative, 
or zero. Though the specifics of the cosmology depend on the initial conditions, there is 
a common nature to the solutions. If dnipi/dt < initially, then the Planck mass starts 
large and quickly drops to some asymptotic value, denoted rhpi. If drUpi/dt > initially, 
then the Planck mass starts out smaller than rhpi and again quickly approaches the 
asymptotic value rhpi. In either CctSG^ ctS long as drUpi/dt is significant, the scale factor 
and the temperature evolve with the changing rupi in a complicated way. Thus, even 
though the energy density in ordinary matter is radiation dominated, the variation in the 
Planck mass alters the dynamics from that of a standard radiation dominated cosmology. 
However, once rUpi veers close to its asymptotic value rhpi, then drUpi/dt ~ and the 
universe evolves in a familiar way. At this point, the descriptions of the cosmology for 
the three possible initial drUpi/dt (positive, negative, or zero) converge. The equations 
of motion reduce to those of an ordinary radiation dominated Einstein cosmology with 
Mo, the usual Planck mass of 10^^ GeV, replaced with rhpi. In particular, this means 
R oc rhpi^^^t^/^, H = l/2t, and dhoriz oc t where R is the scale factor, H is the Hubble 
constant, and cihoriz is the horizon radius. Thus, despite the underlying structure of the 
theory, gravity appears to be described by general relativity with a static gravitational 
constant. 

However, in this early phase of the universe, the gravitational constant can be much 



larger than its value today: rhpi ^ Mq. We call this epoch of large Planck mass the 
MAD epoch. Here we illustrate briefly how a MAD expansion can resolve the horizon 
problem. Once the Planck mass has reached its asymptotic value rhpi, the age of the 
universe scales as t oc rhpi/T'^. A universe which has an early MAD era with large Planck 
mass rhpi » Mo is therefore older at a given temperature than a universe which has 
today's value of the Planck mass Mg for all time. This gives us a hint as to how such a 
scenario may solve the smoothness problem. If the universe is older than in the standard 
model, then much larger regions of spacetime would have come into contact than we had 
previously supposed. This opens the possibility that large regions became smooth without 
violating causality. We will describe this approach to resolving the horizon problem more 
quantitatively later. 

Before we move on to the calculation, we should complete the history of this cos- 
mology. As in a standard cosmology, the energy density in nonrelativistic matter will 
eventually exceed the energy density in radiation. Thus the era of radiation domination 
will end as the universe becomes matter dominated. A matter dominated Brans-Dicke 
cosmology has been well studied^, as have the constraints on such models^. We discuss 
the constraints on MAD expansion with Brans-Dicke gravity in section III.D. We men- 
tion here one of the most severe constraints on our model. As discussed above, in order 
to solve the smoothness problem, our model requires a large value of the Planck mass 
rhpi at some time during the radiation dominated epoch. During the matter dominated 
era, the Planck mass will continue to evolve. However, if the Brans-Dicke parameter is 
greater than 500 as the observations imply, then in pure Brans-Dicke gravity, rUpi will 
not evolve enough during matter domination to reach the value of Mq today. We suggest 
additional physics which may amend this problem, such as a potential for the field or 
a more general scalar theory of gravity. 

We consider a more general class of scalar theories in which the Brans-Dicke param- 
eter u is not constant, still without a potential for the scalar field. We present solutions 



during the radiation dominated era. In this case of variable a;, the Planck mass can 
change more rapidly during the matter-dominated era and could conceivably reach the 
value Mo by today. To check this conjecture, the solutions to the equations of motion 
during a matter dominated era for a general scalar theory would have to be obtained. In 
addition, the constraints on more general models need to be studied in future work. 

Extended inflation^ and hyperextended inflation^ were both developed in the context 
of scalar theories of gravity. In addition to the scalar field which couples to gravity, these 
models require another scalar field, the inflaton field, and a potential for the infiaton. 
The horizon problem is resolved in the usual inflationary way as the scale factor grows 
superluminally during an era of false vacuum domination and then the universe is reheated 
during a period of entropy violation. It is interesting to note that these models also need 
additional mechanism, such as a potential for the Brans-Dicke field, to drive the Planck 
mass down to the value Mo by today. 

In Section II) we present the action and equations of motion for the alternate theories 
of gravity that we are considering. Section III) focuses on Brans-Dicke gravity: III. A) 
presents solutions to the equations of motion during the radiation dominated era, with 
solutions parametrized in terms of the Brans-Dicke field III.B) relates these solutions 
to time evolution; III.C) illustrates the causality condition required to solve the horizon 
problem; and III.D) discusses problems with and constraints on the scenario. Section 
IV) presents a general (albeit somewhat preliminary) discussion of the MAD era in the 
context of more general alternate theories of gravity in which the Brans-Dicke parameter 
u> is not constant. We summarize our conclusions in Section V). 

II) Action 

Brans and Dicke proposed an extension of Einstein gravity in which a scalar field 
usurps the role assumed by the gravitational constant in the Einstein action; that is, in 
Brans-Dicke gravity, the gravitational constant G is not a fundamental constant but is 



instead inversely proportional to a scalar field. More generally, G may be some more 
complicated function of a scalar field ijj: = m^^ = The most general scalar- 

tensor theories^ were originally studied by Bergmann^ and by Wagoner^. Regardless of 
the specific form of the action for such an extension of general relativity is 
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where we have used the metric convention (— , +, +, +), IZ is the scalar curvature, ^matter 
is the Lagrangian density for all the matter fields excluding the field -0) and V{^{^)) is 
the potential for the field if). The parameter u is defined by a; = Stt (^g^Jg^y ■ 

Stationarizing this action in a Robertson- Walker metric gives the equations of motion 
for the scale factor of the universe R{t) and for 

^ „^:ri Stt . „ , dU 1 div 
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where 
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a$ (3 + 2uj) 

U effectively acts as a potential term in the equation of motion for ^. H = R/R is the 
Hubble constant, while p is the energy density and p is the pressure in all fields excluding 
the field. 

The energy-momentum tensor of matter, Tj^^tter; conserved independently of the 
energy-momentum tensor for the scalar, T^^ . The conservation equations are 

^matter :z/ = (5) 



and 



Equation (|^) returns the equation of motion of (J^). It can be shown that in an isotropic 
and homogeneous universe, the /j, — component of equation (|5|) gives 

p = -ip + p)3H . (7) 

Consider the radiation dominated era where p = {tt'^ /30)g^,T'^, p = p/3, and g*{t) is the 
number of relativistic degrees of freedom in equihbrium at time t. Since conservation of 
energy-momentum in ordinary matter does not involve $, we can deduce from eqn (|^) 
that the entropy per comoving volume in ordinary matter, S = {p+p)V/T, is conserved. 
For convenience we define 

s = i^^T^ (8) 

where S = S{4/3){tv'^ /30)gs and gs is the number of relativistic degrees of freedom 
contributing to the entropy. For practical purposes we can take gs = 9*- 

Once the equation of state, and the forms of $(V') and V{'iIj) are specified, these 
equations describe the evolution of the scale factor, the energy density, and the Planck 
mass. In specific, the equations of motion (||) and (J^) and the conservation equation 
determine p{t), and R{t) up to four constants of integration. Notice that S is 

the constant of integration from integrating the energy equation (|^). We take the other 
three constants to be the initial value of $, the constant C defined in eqn. (10) (C oc $), 
and the constant of integration $ given in eqn. (15) (the asymptotic value of $ in the 
radiation dominated era). Given these four initial conditions, the entire cosmology is 
specified: the equations of motion uniquely determine R{t), $(t) and p{t) for all time. 
In contrast, the standard cosmology with a constant Planck mass requires only that the 
value of the Planck mass and two boundary values be specified. The two values needed 
could, for example, be the entropy and the initial value of the scale factor. 

To illustrate how the underlying structure of the Planck mass can alleviate the 
horizon problem, we will present two simple possible scenarios here; that is, we treat 



p{p) = p/3. We first treat the original Brans-Dicke proposal of $ = ^V'^ with oj, 
defined above, constant. The second case we consider is general ^{4>) for which lv is not 
constant. In both cases we take V{iIj)=0, although in another paper we treat the model 
with a nonzero potential. Still, we want to stress that an explanation of the smoothness 
of our observable universe in a cosmology with a variable Planck mass is more general 
than the specific models we study. The crucial ingredients are forms of $('0) and V{ip) 
which resolve the horizon problem without requiring an epoch of vacuum domination nor 
of entropy violation. 

Ill) Case a: a;=constant 

As a first example, take the original Brans-Dicke model where lo is constant and 

mj, = $ = ^-^r . (9) 

Also, we take = 0. In the limit a; = oo, Brans-Dicke gravity is the same as Einstein 

gravity; here we consider arbitrary u> and comment on experimental bounds on u> below. 
We are interested in the behavior of the solutions during the hot radiation dominated era 
of the early universe. We assume the matter energy density is negligible. Then p = p/3 
and p — 3p=0. 

We first obtain solutions to the equations of motion. Instead of finding R{t), 
and T{t), it is more tractable to parameterize R, T, and hence H and the horizon radius, 
c^horiz, by We then find approximate solutions for $ as a function of time in different 
regimes. Before moving on to the second case of a; 7^ constant, we address the horizon 
problem in the context of our solutions. 

III. A) Solutions to the Equations of Motion 

With (jj constant and V{'ip) = 0, the $ equation of motion during the radiation 
dominated era reduces to $ -|- 3H^ = 0, so that 

^R^ = -C , and iy = --S- , (10) 



where C is a constant of integration which can be positive, negative, or zero. If C = 0, 
then $ = 0, the Planck mass has a constant value which we may call rhpi, and the 
universe evolves in the usual radiation dominated fashion, but with G oc l/rhpi. Note 
that if C > 0, then $ < 0, while if C < 0, then $ > 0. Here we take k = to illustrate 
the behavior of the solutions. [In the appendix, we present the solutions for n — ±1.] 
First we solve equation for H: 
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where 7(t) = (87r/3)(7r^/30)(74<(t). Note that all three terms inside the square root are 
positive quantities. We choose the sign in front of the square root in such a way as 
to obtain an expanding universe with H > 0. Thus, for C > we take the + sign in 
equation (JTj), whereas for C < 0, we take the — sign in equation (JT]). Throughout the 
rest of the paper, the upper sign in equations will refer to the case C > and the lower 
sign to the case C < 0. Substituting $ from eqn. (^Ol) into the square root in eqn. ( [TT|) 
and using H = R/R, we have 

1/2- 
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We define 



= 4S^/^C- (^^) (13) 
and note that x is always a real positive quantity. The integral of eqn. (|12D becomes 

dx' 
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where subscript i refers to initial values. We find the (positive) solution 
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Here, $ = $i exp 



-(l/e)arcsinh(x5 
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i.e. we have absorbed the constants of inte- 



±(l-|-2a;/3) 



1/2 



(as noted above, an expanding universe corresponds to the + in e if C > 0, or the — in e 
if C < 0). For convenience, we define 



e = eln 



(16) 



In many of the equations below we express the field in terms of rather than $. As 
we will show below, G is always positive semidefinite for any value of C: it ranges from 



= for $ = $ to O = cxD for $ far from $. From eqns. (|T3D , (^5]) and (|T^) , we obtain 
an expression for the scale factor. 
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Note that the product (Ce) is always positive semidefinite. 



It follows from adiabaticity that 
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The Hubble constant -ff($) is obtained from 



where in the last step we used eqn. (^Up . We find 

sinh^ G {sinh O + 2e cosh 0} . 
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The comoving horizon size is 
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(18) 



(19) 



(20) 
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where we have used $ = —C/R{^)^ in the last equality. We can integrate this to find 



If $(t = 0) starts out far from $, i.e. G(t = 0) ^ 1, eqn. (|22|) becomes 

(ihoriz f C \ e exp [-6] 



(23) 



As discussed above, -R($), and (ihoriz($) are determined only up to 

the arbitrary constants C, and 5^. For instance, eq. (^Sp shows that one can choose 
the temperature at a given value of $ by choosing C /S and $ appropriately. The fourth 
and last constant of integration, $(t = 0), is determined when is found in section 
III.B. 

Even before we determine we can understand the general sketch of the uni- 

verse's evolution. We will find that, in all cases, the field $ asymptotically approaches 
the value $: for C > 0, $ approaches $ from above; whereas for C < 0, $ approaches $ 
from below. [For C = 0, $ = $ for all time]. 

Let us first consider the case of C > 0. As mentioned previously, this corresponds to 
$ < and e > 0. In order for the scale factor to satisfy -R > 0, from (^7]) we can see that 
we need O > 0, i.e., $ > $. In addition, we need the scale factor to grow in time; again, 
this requires $ < 0. In short, for C > 0, $ starts larger than $ and decreases towards 
the asymptotic value $. 



For the case of C < 0, we have $ > from (^0|) and e < 0. To obtain i? > 0, we need 
> 0, which in this case corresponds to the opposite limit of $ < One can show that 
as $ grows towards its asymptotic value (i.e. drops), as long as |e| > 1/2 (i.e. u > 0), 
dR/dO < 0; the scale factor grows in time. In short, for C < 0, $ starts smaller than $ 
and grows towards the asymptotic value $. 

As we have seen, as $ approaches $, for |e| > 1/2 (w > 0), -R($) grows and thus the 
temperature drops adiabatically. In addition, one can show (again for |e| > 1/2) that the 
comoving horizon size grows, as does H~^R~^. We'll see below that the size of a causally 



cools below the temperature of matter radiation equality (Teg ~ 5.50^^^ eV, where Qm 
is the fraction of the critical density contributed by matter), the reign of radiation yields 
to that of matter and the nature of the solutions changes. During matter domination, 
the equation of state is p{p) = 0, and p — 3p = p. This alters the dynamics considerably. 
Thus there is a built in off-switch to end the radiation dominated behavior of -R($), T($), 
and 4oriz($)- 

We will quantify these statements below and find constraints on some of the con- 
stants of integration needed to resolve the horizon problem. In the next section, we find 
approximate descriptions of the behavior of $ as a function of t. Actually, the resolution 
of the smoothness problem and the evolution of the cosmology can be understood with- 
out knowing $ as a function of t. The universe will pass through the familiar stages of 
baryogenesis, nucleosynthesis, matter domination etc. as the temperature passes through 
the relevant energy scale. To follow the evolution of the universe all one needs to know 
is the temperature as a function of $. One need not know the actual age of the universe. 
Still, to ground the solutions in a slightly more familiar setting we will indicate below 
how the universe evolves in time. 

III.B) The Age of the Universe 

We will determine the time evolution of the Brans-Dicke field $ in two different 
limits: $ far from $ (G ^ 1), and $ ~ $ (0 ^ 1). As we have seen, initially $ may be 
large for $ < (C > 0) or small for $ > (C > 0). While $ is far from its asymptotic 
value $, the term $/$ contributes significantly to the equations of motion (see eqn. (^) ), 
and the evolution of the universe is modified relative to that of an Einstein universe in a 
complicated way as eqns. (^7[)-(|2^) show. Once $ ^ $ however, the universe evolves with 
time as an ordinary radiation dominated cosmology with the Planck mass Mo replaced 
with rhpi = 

To uncover $(t), return to $ = —C/R^ (c.f. eqn. (p!0| ) ). We integrate this equation 



to find 

/ d^'R^i^') = -Ct , (24) 

where -R($) is given in eqn. (^7]) . To get a rough feehng for how $ changes with t we 
find approximate solutions to this integral for two regimes: (1) $ far from the asymptotic 
value $ and (2) $ ^ 

(1) $ far from $ 

First we consider the early regime where $ is far from $, i.e. 3> 1. The integral on 
the left hand side of eqn. ( p^ is easiest to evaluate if rewritten in terms of O rather than 
$. Then we can approximate sinhG ~ e®/2 in evaluating the integral. The lower limit 
of the integral is determined by the boundary condition: R(t = 0) ^ 0. For all values 
of C, this initial value of the scale factor requires 0(t = 0) ^ oo; i.e., $(t = 0) ^ oo for 
C > while $(t = 0) ^ for C < 0. Thus, our fourth and last integration constant is 
determined. [Given C, S, and $, R{t = 0) and $(t = 0) contain the same information via 
eqn. ( pTTD ]. We can now evaluate eqn. ( pif ) to find 
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We see from eqns. ( [25| ) and ( p!7D that initially R{t) oc Remember 
e = ±1/2(1 + 2a;/3)^/^. It is interesting to consider the nature of these solutions for 
large deviations from Einstein gravity (i.e. small uj). Take C < for a; ^ 0, and thus 
e —1/2. In this limit, ^ —>■ t while R —>■ constant. Note that this behavior of the scale 
factor could also be seen directly from eq. (|l^ . On the other hand, for C > with 
^ 0, e ^ +1/2 and $ f"^/^ while R t^/^. Again, the behavior R oc could 
be seen directly from eqns. (|T7| ) and (ITSl). 

(2) $ ^ $ 

The previous approximation breaks down for $ $. Again, it is easiest to work 



order contribution to the integral yields 



or, equivalent ly, 
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(26) 



(27) 



2^2^3/2$l/2^^ 

where the plus sign refers to C > and the minus sign to C < 0. [By assumption, we 
are working near $ $ so that the exponent must be small for this approximation to 
be valid.] 

Since eqn. (^7|) implies that 0(t) oc t~^l'^ (times a positive constant), we have 
i?(t) oc oc go^ as $ approaches $, the universe evolves as an ordinary radiation 
dominated universe with one modification; the Planck mass Mo is replaced by $^/2 

In the standard Hot Big Bang model described by Einstein gravity, the age of the 
universe as a function of temperature is given by teinst = 2^f/2r2 ■ As $ approaches $, 
we can see from eqns. (|27|) and (|T8| ) that the age of the universe as a function of T($) 
mimics this form, 



m 



rripi 



(28) 
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where nipi = $^/2 Incidentally, this is exactly the result one obtains for C = ($ = 0) 
and $ = $. As discussed in the introduction, the universe is older for a given temperature 
if rhpi > Mo than an Einstein universe with Planck mass Mq. We will show below that if 
the comoving horizon volume is to become smooth at high temperatures, then rhpi must 
greatly exceed Mo- 

III.C) Horizon Condition and Discussion 

We show here how a Brans-Dicke, radiation dominated cosmology can resolve the 
horizon problem if there is an early MAD epoch with rhpi ^ Mo- 



To explain the smoothness of our present universe, a region causally connected at 



Since we can only see back to the time of decoupling, or perhaps nucleosynthesis, we 
can take the present comoving Hubble radius, l/{HoRo), as a measure of the comoving 
radius of the observable universe; here Ho is the Hubble constant today and Ro is the 
scale factor today. Then the smoothness of the observable universe can be explained if a 
comoving region of radius at least as large as 1/HoRo is in causal contact at some time 
tc before nucleosynthesis, i.e., 

>ir-Er ' (29) 



HcRc RqHo 

where subscript c denotes values at the time causality is satisfied. 

We can express both -ff ($) and Ho in terms of the temperature and the Brans-Dicke 
field $. Substitution of expression (^8D for T($) into equation (|0D for gives 



=7V2^^1|sinhe + 2ecoshe} . (30) 



The Hubble constant today can be written as 



^2 



H.-al'^jl (31) 

where To = 2.6 x 10"^"^ GeV, Mo is the value of the Planck mass today, and cto = 
7{to)'rjo = (8vr/3)(7r^/30)5rH<(^o)^o5 where rjo ~ 10^ — 10^ is the ratio today of the energy 
density in matter to that in radiation. Also, we use adiabaticity, RT = S^^^ oc [S/ g^Y^^, 
to write the causality condition as 

2e ^ ^Mo 

Tc sinhGc + 2ecoshOc ~ To 

where 13 = {'^{tc) / GioY^'^{g*{tc) / g*{to))~^^^ ■ To resolve the horizon problem, this con- 
straint must be satisfied prior to matter /radiation equality. 

Although it is possible for the causality condition (^2]) to be satisfied while $ is still 
far from $, we find that the solution to the horizon problem that deviates the least from 



1/2 ~ 

value of $c oc mpi{tc) that solves causality is given by $c — ^- [For C < 0, $ is always 
less than $, yet the previous statement still holds to better than 1 %; if $ <^ $, both $ 
and $ are driven to higher values than if they are equal to each other] . From now on, we 
will examine the causality constraint for $ near its asymptotic value $. 

For $ ~ $, ~ 0, sinhG ~ 0, cosh© ~ 1, and the causality condition becomes 
simply 

^ > 3^ (33) 

where mpi{tc) ~ rhpi = We can specify the temperature at which we would like to 

resolve the causality dilemma. We are free to choose the temperature at which $ = $ 
since this is equivalent to making an appropriate choice for the ratio of the arbitrary 
constants S/C (see eq. ([T8|) ). [Since S/C oc T"^/$, this amounts to making a choice for 
$(tc); in principle one should check that this choice is consistent with measurements of 
G/G today. However, since $ oc R~^, in many cases the time derivative may be quite 
small and therefore unobservable by the present epoch.] 

As an example, we consider T($c ~ $) ~ 1 MeV, roughly the temperature of 
primordial nucleosynthesis. Then condition (|33|) requires 

mpi> 10'' Mo . (34) 

We can verify that the smoothness problem is explained by an old universe. We 
showed with our approximate expressions for nipi as a function of time, that when $ $, 
the universe evolves as an ordinary radiation dominated universe with one modification; 
Mo is replaced by rhpi. In this limit 

t($) = !!h±^ ^35) 

27V2T($)2 ^ ^ 

Since tc oc I /He and to oc I /Ho, eqn. ( p9D is equivalent to the statement that 

^ -B- ■ (36) 



Since rhpi does in fact exceed Mq we see from eq. (^5]) that the universe is older at 
a given temperature than in the standard cosmology. Writing eq. (|35| ) in terms of teinst, 
the age of a cosmology described by Einstein gravity, gives 

m = (^) (37) 

at a given temperature. In a standard cosmology, teinst ~ (MeV/T)^ sec. So, at 1 MeV, 
^einst ~ 1 sec and ~ 10^ sec ~ 3 yrs. 

If, instead, we take T($c ~ to be the temperature of matter /radiation equality, 
about 5.5 VtMh"^ eV, then the causality condition requires 

mpz>O(10')M„ . (38) 

At T ~ 1 eV, teinst ~ 10^2 sec ~ 10^ yrs, and t($) ~ 10^ teinst ~ 10^ yrs. We see that 
the smoothness of the observable universe is resolved in this MAD model by aging the 
universe. 

III. D) Problems and Constraints 

The obvious difficulty with this resolution to the horizon problem is fixing the value 
of the Planck mass to be Mo by today. In the Brans-Dicke model studied here without a 
potential, the Planck mass will be hard pressed to make it to the value Mq today. During 
the matter dominated era, $ will initially continue to decrease with time for C > 
and increase with time for C < 0.^ For C < then, the Planck mass will only grow 
larger. For C > 0, it is conceivable that $ will approach the value of during the 
matter dominated era. However, observations constrain the parameter w to be > 500 for 
a massless Brans-Dicke theory. The rate at which $ changes depends on uj and is very 
suppressed for large w. Thus a large to would confine $ to near its value at the time 
of matter radiation equality, which, as we have seen, may be large. For example, with 
UJ = 500 and ^^/^ = lO^Mo at ~ 1 eV, then today > 80Mo. This hmit can 

2 



interactions measured in the time-delay experiments fall off rapidly outside the range 
over which the $ field acts. If $ has an associated mass m^, then the range over which 
$ acts A ~ could be smaller than the distances over which the observations are 

sensitive. Therefore a massive Brans-Dicke model would elude observation even if cv is 
small^. However, in another paper ^° we discuss the case of Brans-Dicke with a potential 
and find that it is diflacult to simultaneously satisfy bounds on time-delay experiments 
as well as other constraints. 

There are also observations of the rate of change of the gravitational constant. These 
observations impose a much weaker constraint than the time delay experiments. They 
suggest u > b. For comparison, if a; = 5 and $-'^/^ = IO^Mq, then one could have 

^ Mo up to 0(1). 

Another issue of concern is the value of the Planck mass, and thus the Hubble 
constant, during nucleosynthesis. If nipi ^ Mg during nucleosynthesis, then the predicted 
elemental abundances will be affected. To resolve the horizon problem, we found that 
the asymptotic value rhpi = had to be much larger than Mq. If, for example, 

rupi ~ rhpi » Mo during nucleosynthesis, then H{^) oc ^_ [f^ and the large Planck mass 
slows the expansion of the universe. Consequently, the temperature at which the weak 
interactions freese out is lowered, the n/p ratio is maintained at its equilibrium value 
longer, and the value of the n/p ratio during nucleosynthesis is smaller. This works 
to decrease the production of He^. Compatibility with observations would then force 
Qb, the fraction of critical density in baryonic matter, to be larger. Actually, since the 
Hubble constant can be so much smaller than in the standard model, a situation may 
arise where the weak interactions are still in equilibrium during the time of He^ synthesis; 
then the nucleosynthesis calculations would have to be redone. We have not investigated 
the consequences for abundances of other elements (such as deuterium, lithium, and He^). 
We suspect that, unless the Planck mass has returned to its present value by the time of 
nucleosvnthesis. matching observations on all elements simultaneouslv will be impossible. 



One could insist that the causahty condition is solved for temperatures greater than 
an MeV and then invoke a potential to drive nipi to Mq by the time of nucleosynthesis. 
This would also accomodate the C = scenario where $ is constant at the value rhpi 
needed to solve causality. For the previous results to hold, the potential would have 
to remain inconspicuous during the early evolution. Again, we treat the model which 
includes a potential in a separate paper and find that it is difficult to satisfy all the 
constraints. 

Next we will consider a more general scalar type theory where $('(/') is such that 
is not constant. We treat this case below. 

IV) Case b: u T^constant: 
IV. A) Solution to the Equations of Motion 

Here we extend the analysis to the more general case of w not constant, again with 
no potential for the Brans-Dicke field. We will find the solutions to and 
T($) here. The solutions during the radiation dominated era are very similar in spirit to 
the previous solutions for uj = constant. However, the case of w not constant does allow 
the possibility of a small value of w at early times that matches onto lo > 500 today. 
Although we have not yet worked out the evolution for specific examples of changing uj 
during the matter dominated era, one can hope that $ can change rapidly enough during 
the period of matter domination to reach the value today. For certain choices of $(^/'), 
uj{i(j) will grow in order to match the observational constraints of a; > 500 today without 
invoking a potential. (Actually, the observational constraints need to be reinterpreted if 
u 7^ const ant.) 

If $ has any functional form other than the minimal oc "0^, then uj = 8tc (^g^^g^y will 
also be a function of t/j. The equation of motion (Q) in the radiation dominated era for 
the case of no potential is 



so that 

^R' = -7 ^VttT^ and H = -X - , ^^^^ , , . (40) 

(l + 2a;/3)i/2 3$ 6(1 + 2a;/3) ^ ^ 

We can use results (^Of) and (^) in eqn (^) as we did in the previous section. This time, 
we define 

X = 4S^/^'yC-^^R^ (41) 

(this time, there is no factor of (1 + 2(^/3)"^ in our definition of x). Again, we obtain 



eqn. (Ill) . Since a;($) is as yet unspecified here, we define the right hand side of eqn. 



(0) to be 

E2($) = ± /^^ ^(1 + 2^73)^/' . (42) 
Then the solution to (|T^ ) in this case becomes 

X"^/^ = sinh[E2($) + arcsinhx,"^^^] . (43) 

From eqns. (^l]) and (^) , we can find 

= ^ 1 (44) 

2^2/3^1/2 $i/2ginh[E2($) + arcsinhx,"^/^] 

— 1/2 

From now on we define cr^ = arcsinhXj and absorb this term into the definition of 
E($) = E2($) + a,. 

It follows that 

T($) = ^^?I^$i/2ginhE . (45) 
For use in the constraint (^^j), we find H{^) as before, 

iy($)=47^/^(§)'(ff^^^^ • (46) 

iy($) can be expressed in terms of the temperature and E using eq. (|45| ) in eq. (^6|). 

. .\,...J sinhE + (l + 2a;/3)V2coshE} . (47) 



The comoving horizon size is 



1 1 



tanh E tanh a, 



(48) 



In order to solve the smoothness problem, we need )r{^ ) ~ h^r before. 
This gives a constraint similar to the uj =constant scenario. Written in terms of the 
temperature, the constraint is 



T($ J {sinh + (1 + 2ojj2,) V2 cosh T.^} ^ ' 

where (3 is defined below eq. (^2]) and again the subscript c indicates the values at the 
time causality is solved. Again, the causality condition is satisfied with the least deviation 
from Einstein gravity for small (i.e. sinhE^ ~ and coshEc ~ 1); then eqn. (^) 
becomes 

We see that a large early value of the Planck mass is again required to solve the causality 
condition at high temperatures. 

Again, as in the case of pure Brans-Dicke gravity, $ approaches an asymptotic value 
$ from either above or below. We can see this by looking at eqn. (^Dp : for > 0, 
we can see that $<C/i?"^— >Oasi?— >oo, and $ approaches some constant value 
$. For consistency, we want the scale factor in eqn. (|4^) to grow very large as $ 
approaches $. Thus, as $ we can see from eqn. ( p4D that we need E — > 0, i.e., 

f§ (1 + 2a;/3)^'^^ = =Farcsinh(x^ ^/^), where the upper (lower) sign refers to C > (<)0. 
This last relation tells us the asymptotic value $ as a function of initial values Xi c^nd $ j . 

As an example we take $('0) = $exp (^tp/ipj where ip is a constant mass scale and 
$ is a constant with units of mass squared. Then 



The description of this cosmology is very similar to the oj =constant scenario. Ini- 
tially, $ is either positive, negative, or zero. If $ is positive, then the Planck mass starts 
small and increases toward the asymptotic value fhpi = If $ is negative, then the 

Planck mass begins large and decreases toward fhpi. To satisfy the causality condition 
when the Planck mass is at the asymptotic value rhpi, we need ifipi/Mo <; (3T /Tg. 

Although this yields a very similar description of the early universe as the uj constant 
scenario did, there is one advantage. Since u> oc u) is small for large values of For 
small a;, $ changes rapidly with time when the universe is dominated by nonrelativistic 
matter. This buys $ the chance to move toward by today. As $ decreases, uj increases 
and automatically turns off the change in $. A more thorough treatment of the equations 
of motion during matter domination with u) T^constant is needed to verify this suggestion. 

V) Conclusions 

During a MAD era, the Planck mass is large and the universe is older at a given tem- 
perature than in a standard cosmology. The comoving size of a causally connected region 
1/HR is correspondingly larger. Since the observed universe only reaches out to about 
recombination (or possibly back to nucleosynthesis) , a rough estimate of the comoving ra- 
dius of the observable universe is the present comoving Hubble radius H~^/Ro ~ Ato/Ro, 
where Atg is roughly the time elapsed since the universe became matter dominated. The 
smoothness of the observable universe can be explained if the universe ages sufficiently 
so that the comoving horizon size at the time of nucleosynthesis is H~^/R > Ato/Ro- 
Observations of the age of the universe from the Hubble diagram, nucleocosmochronol- 
ogy, and ages of globular clusters only place limits on the age of the universe subsequent 
to the time when stars formed; thus, the universe may in fact become much older in the 
radiation dominated era than one would expect from the standard model. 

We found that a MAD era with large Planck mass ages the universe sufficiently 
to resolve the horizon problem. For a Brans-Dicke theory without a potential during 



the radiation dominated era, the Planck mass approaches an asymptotic value, rhpi. 
This asymptotic value can be chosen to satisfy the causality condition at a specified 
temperature T^. rhpi/Mo > TcJTq. For instance, for causality to be satisfied by the 
time of nucleosynthesis, we need the asymptotic value rhpi > lO^Mc,. Or, for causality 
to be satisfied by the time of matter/radiation equality, we need ihpi > lO^Mo. If the 
Brans-Dicke parameter uj is greater than 500 and constant, rupi will be unable to reach 
the value Mo = 10-^^ GeV by today. Also, if the Planck mass is large during primordial 
nucleosynthesis, then the standard big bang nucleosynthesis predictions will be altered. 
We suggest that rUpi may be driven down to almost its present value by the time of 
nucleosynthesis if there is a potential in the theory for $ (however, in a separate paper^° 
we show that such an approach is typically overconstrained) or if u) is allowed to vary 
as a function of time. (In a separate paper we also discuss the fiatness problem and the 
monopole problem. 

For a more general scalar theory of gravity than the Brans-Dicke proposal, where u 
is not constant, we find again that the Planck mass can approach an asymptotic value as 
the universe evolves. As before, if rhpi /Mo ^ Tc/To, then the smoothness of the present 
observable Hubble volume can be explained in this model. In addition, if uj ^^^constant, 
then the Planck mass may change rapidly enough with time to reach the value of Mo 
by today. The equations of motion during the matter dominated era for a more general 
scalar theory need to be studied in detail, as do the constraints on general models. 
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Appendix 



We present here the solutions to the equations of motion during the radiation dom- 
inated era for a Brans-Dicke theory with k = ±1. The $ equation of motion reduces to 
$ + = so that, as before, 



-C and H 



3$ 



(52) 



Solving equation for H with C 7^ and k 7^ gives 



H 



2$ 



l + 2eVl + X-Q2x2 



(53) 



where e = '^^^ — , x is defined as in (p!3| ) and 



7 



2^8/3 



K 



(54) 



Using if = R/R, the definition of x, and rearranging, we are left with the integral 

dx' 



/2 



(55) 



3 y $' 

Integrating this equation and using R = (eC/52/3) (7$) -1/2^1/2^ from the definition of 
X, we find. 



R 



eC 1 



1 



1/2 



(56) 



^2/3^1/2 $1/2 1^ gj^j^2 Q^Q2 exp (-26) 

where, as before, O = eln($/$). The temperature of the universe is found from adia- 
baticity to be 

^ (57) 



The causality condition becomes 



1/2 



-2e 



(sinh^ + exp(-2ec))^/2 



sinh^ + 2e sinh 6^ cosh 6^ + Q^{1- 2e) exp(-2e. 



rp 
J- n 



(58) 



Notice that as (pq), (pTl), and (^8[) reduce to the corresponding results for a flat 



in the case of the flat universe. We wiU further discuss the case of nonzero curvature, 
together with the issue of flatness, in another paper. 
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